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Abstract 

The procedure for obtaining integrable open spin chain Hamiltonians via reflection matrices is 
explicitly carried out for some three-state vertex models. We have considered the 19-vertex models 
of Zamolodchikov-Fateev and Izergin-Korepin, and the Z2-graded 19-vertex models with s/(2|l) and 
osp{l\2) invariances. In each case the eigenspectrum is determined by application of the coordinate 
Bethe Ansatz. 
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1 Introduction 

One-dimensional quantum spin chain Hamiltonians and classical statistical systems in two spatial dimen- 
sions on a lattice (vertex models), share a common mathematical structure responsible by our under- 
standing of these integrable models [1, 2, 3]. If the Boltzmann weights underlying the vertex models are 
obtained from solutions of the Yang-Baxter (YB) equation the commutativity of the associated transfer 
matrices immediately follow, leading to their integrability. 

The Bethe Ansatz (BA) is a powerful method in the analysis of integrable quantum models. There are 
several versions: coordinate BA [4], algebraic BA [5], analytical BA [6], etc. The simplest version is the 
coordinate BA. In this framework one can obtain the eigenfunctions and the spectrum of the Hamiltonian 
from its eigenvalue problem. It is really simple and clear for the two-state models like the six-vertex 
models but becomes tricky for models with a higher number of states. 

The algebraic BA, also known as Quantum Inverse Scattering method, is an elegant and important 
generalization of the coordinate BA. It is based on the idea of constructing eigenfunctions of the Hamil- 
tonian via creation and annihilation operators acting on a reference state. Here one uses the fact that 
the YB equation can be recast in the form of commutation relations for the matrix elements of the mon- 
odromy matrix which play the role of creation and annihilation operators. From this monodromy matrix 
we get the transfer matrix which commutes with the Hamiltonian. 

Imposing appropriate boundary conditions the BA method leads to a system of equations, the Bethe 
equations, which are useful in the thermodynamic limit. The energy of the ground state and its excita- 
tions, velocity of sound, etc., may be calculated in this limit. Moreover, in recent years we witnessed an- 
other very fruitful connection between the BA method and conformal field theory. Using the algebraic BA, 
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Korepin [7] found various representations of correlators in integrable models and more recently Babujian 
and Flume [8] developed a method from the Algebraic BA which reveals a link to the Gaudin model, ren- 
dering solutions of the Knizhnik-Zamolodchikov equations for the SU{2) Wess-Zumino-Novikov-Witten 
conformal theory in the quasiclassical limit. 

Integrable quantum systems containing Fermi fields have been attracting increasing interest due to 
their potential applications in condensed matter physics. The prototypical examples of such systems arc 
the supersymmetric generalizations of the Hubbard and t-J models [9]. They lead to a generalization of 
the YB equation associated with the introduction of the a Z2 grading [10] which leads to appearance of 
additional signs in the YB equation. 

When considering systems on a finite interval with independent boundary conditions at each end, 
we have to introduce reflection matrices to describe such boundary conditions. Integrable models with 
boimdarics can be constructed out of a pair of reflection JC-matrices K^{u) in addition to the solution 
of the YB equation. Here K~{u) and Ar+(u) describe the effects of the presence of boundaries at the left 
and the right ends, respectively. 

Integrability of open chains in the framework of the quantum inverse scattering method was pioneered 
by Sklyanin relying on previous results of Cherednik [1 1] . In reference [12] , Sklyanin has used his formalism 
to solve, via algebraic BA, the open spin- 1/2 chain with diagonal boundary terms. This model had already 
been solved via coordinate BA by Alcaraz et al [13]. The Sklyanin original formalism was extended to 
more general systems by Mezincescu and Nepomechie in [14]. 

In this paper we consider the coordinate version of the BA for the trigonometric three-state vertex 
models with a class of boundary terms derived from diagonal reflection ii'-matrices. These models are 
well-known in the literature: the Zamolodchikov-Fateev (ZF) model or ^l^^'' model [15], the Izergin- 
Korepin (IK) model or A\ model [16] and two .Z2-graded models, named the sZ(2|l) model and the 
osp(l|2) model [17]. 

In the context of the coordinate BA, we propose here a new parametrization of wavefunctions. This 
result is important since it allows us to treat these 19-vertex models in the same way. Moreover the 
coordinate BA for these threes-states models becomes simple as for two-state models terms [13]. 

The main goal in this paper is to reveal the common structure of these 19-vertex models with boundary 
terms which permits us to apply the BA method, unifying old and new results. 

The paper is organized as follows: We introduce the algebraic tools in Section 2. In section 3, we apply 
the coordinate BA method for a general open chain Hamitonian associated with four 19-vertex models. 
In sections 4,5,6 and 7 the energy eigenspectra and the corresponding Bethe equations are presented for 
each model. In section 8 we discuss about the graded and non graded solutions for 19-vertex models. 
Section 9 is reserved for the conclusion. 

2 Description of the model 

To determine an integrable vertex model on a lattice it is first necessary that the bulk vertex weights be 
specified by an 7?.-matrix TZ{u), where u is the spectral parameter. It acts on the tensor product (g) V'^ 
for a given vector space V and satisfy a special system of functional equations, the YB equation 

ni2{u)ni3{u + V)n23iv) = T123{v)TIi3{u + V)nx2{u), (2.1) 

in «) y2 ^ V^, where 7^l2 = 7^ (g) 1, 7^23 = 1 «) 7^, etc. 

An TZ matrix is said to be regular if it satisfies the property 7^(0) = P, where P is the permutation 
matrix in ® V^: P{\a) ® ]/?)) = ]/3) ® \a) for \a) , e V. In addition, we will require [14] that TZ{u) 
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satisfies the following properties 



regularity 
unitarity 
PT — symmetry 
crossing — symmetry 



7^l2(o) = /(o)l/2Pl2, 

Tli2{u) = U^n%{-u - p)Ur\ (2.2) 



where f{u) = xi{u)xi{—u), U denotes transposition in the space i , p is the crossing parameter and U 
determines the crossing matrix 

M = U*U = M*. (2.3) 

Note that unitarity and crossing-symmetry together imply the useful relation 

MiTelK-u - p)Mf i7e^2(« -p) = f{u). (2.4) 

The boundary weights then follow from /C-matrices which satisfy boundary versions of the YB equation 
[12, 14]: the reflection equation 

■Ri2{u - v)K^{u)n{\*^{u + v)K^{v) = K^{v)Hi2{u + v)K^ {u)Tl{''^^ {u - v), (2.5) 

and the dual reflection equation 

Tli2{-u + v){K+Y'{u)M^^n{\'%-u -V- 2p)Mi{K+Y^{v) 

= {K+Y^{v)Mi'Ri2{~u~v-2p)M-\K+f'{u)n'^^^{-u + v). (2.6) 

In this case there is an isomorphism between K~ and K'^ : 

K-{u) :^ K+{u) = K-{-u- pfM. (2.7) 

Therefore, given a solution to the reflection equation (2.5) we can also find a solution to the dual reflection 
equation (2.6). 

In the framework of the quantum inverse scattering method, we define the Lax operator from the TZ- 

matrix as Laq{u) = TZaq{u), where the subscript a represents auxiliary space, and q represents quantum 
space. The row-to-row monodromy matrix T{u) is defined as a matrix product over the N operators on 
all sites of the lattice, 

T{u) = LaN{u)LaN-liu) ■ ■ ■ ial(w). (2.8) 

The main result is the following: if the boundary equations are satisfied, then the Sklyanin's transfer 
matrix 

t{u) = Tra iK+{u)T{u)K-{u)T-\-u)) , (2.9) 

forms a commuting family 

[tiu),t{v)]=0, yu,v (2.10) 

The commutativity of t{u) can be proved by using the unitarity and crossing-unitarity relations, the 
reflection equation and the dual reflection equation. It implies integrability of an open quantum spin 
chain whose Hamiltonian (with K~{0) = 1), can be obtained as 



' ^ , troi^+(0)g^,o 



2 du 
fe=i 



M=0 



trii'+(0) 
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and whose two-site terms are given by 



(2.12) 



■u=0 



in the standard fashion. 

Here we will extend our discussions to include the Z2-graded vertex models. Therefore, let us describe 
some useful informations about the graded formulation. 

Let V = Vb® Vi be a Z2-graded vector space where and 1 denote the even and odd parts respectively. 

s 

Multiplication rules in the graded tensor product space V V differ from the ordinary ones by the 

s s 

appearance of additional signs. The components of a linear operator A i^i B £ V V result in matrix 
elements of the form 

(A ^ B)l'p = (-)P(/3)(P(«)+P(7)) A^^Bps- (2.13) 

S S 

The action of the graded permutation operator P on the vector \a) ^ \P) G V <^ V is defined by 

V \a) ^ 1/3) = (-)P(«M/3) 1/3) ^ \a) =^ [VfJ^ = (-)f(«)f(^)<5„5 5^^. (2.14) 
The graded transposition st and the graded trace str are defined by 

(^'*)a/3 = (-)(^("'+^)^(^)A;3„, str^ = ^(-)f(«)A„„. (2.15) 

a 

where p{a) = 1 (0) if \a) is an odd (even) element. 

For the graded case the YB equation and the reflection equation remain the same as above. We only 
need to change the usual tensor product to the graded tensor product. 

In general, the dual reflection equation which depends on the unitarity and cross-unitarity relations 
of the 7?.-matrix takes different forms for different models. For the models considered in this paper, we 
write the graded dual reflection equation in the following form [18]: 



= {K+y*^{v)Mrnf^'''{-u -V- 2p)M^\K+y*' {u)ni\''*'{-u + v), 



(2.16) 



and we will choose a common parity assignment: p{l) = p{3) = and p{2) = 1 , the BFB grading. 
Now, using the relations 



TZf2''*'{u) = hR2i{u)h, 7efi^'*=(tt) = /iiii2(u)/i and IK+{u)I = K+{u) 



(2.17) 



with I = diag(l, —1, 1) and the property [MiM2,TZ{u)] = we can see that the isomorphism (2.7) holds 
with the BFB grading. 

The three-state vertex models that we will consider are the Zamolodchikov-Fateev (ZF) model, the 
Izergin-Korepin (IK) model, the sZ(2|l) model and the osp{l\2) model. Their 7?.-matrices have a common 
form 



n{u) = 



X2 



X3 




2/7 



X5 



ye 



\ 



xr 



X3 



X2 



(2.18) 



Xl I 
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satisfying the properties (2.1-2.4) together with their graded version. 

In order to derive the bulk Hamiltonian, it is convenient to expand the normalized i?-matrix {R = PTZ) 
around the regular point u = 



R{u, r/) = 1 + u{a-'^H + pi) + o(u^), 



(2.19) 



with a and f3 being scalar functions. Therefore Hk,k+i in (2.11) is the H in (2.19) acting on the quantum 
spaces at sites k and fc + 1 . 

Using a spin language, this is a spin 1 Hamiltonian. In the basis where is diagonal with eigenvectors 
\+,k) ,\0,k) ,\ — ,k) and eigenvalues 1,0,-1, respectively, the bulk Hamiltonian density acting on two 
neighboring sites is given by 





( 


Z5 


Z7 


1 


Z6 




23 






1 




Z5 










Hk,k+i = 






eze 




ezi 




















Z5 




1 








Z3 




Z6 


1 


Zl 


Zh 

Zl 



\ 



(2.20) 



which can be easily written in terms of the usual spin-1 operators: 



Hk,k+i = ez^ + ^{h-Z5)[S^,-SI^,] + ^{z, + -z,-'2ez^)[iS',f + iSUif] 

+ i(22i -Z7- Z7)S^S'k+i + \{^zx +Z7 + Z7 + Aezi - Az^ - A-z^)){SlS'^+^f 

1 

^2;\2 Qz 



k+l) 



'^'^[^k ^k^k+l^k+l'^^k ^k^k+l^k+l'^^k^k ^k+l^k+l+^k^k ^k+l^k+l]- (2-21) 



Here we have used the sign e = ±1 to explicitly include the graded (e = — 1) models for which the tensor 
products in (2.21) are BFB graded and 



10 
5^ = I 

0-1 



1 
5+ = \/2 ( 1 





\/2 ( 1 

1 



(2.22) 



Finally, we can consider the boundary terms of H which are derived from the diagonal solutions = 
dia,g{kf-^^,k22, kf^) of the reflection equations. In general, at the flrst site (left) it has the form 



1 dK-(u) 

; 

2 du 



11 



'22 



'33 



" ~ 2 du 



1,2,3. 



(2.23) 



5 



and for the last site (right) it has the form 




(2.24) 



To compute the term troX(|(0)iJAr,o we have to use (2.21) in order to get .ffjv.o- In practice, it is equivalent 
to take the trace of K^{0)H2i where 

hii hi2 hi3 

(2.25) 

/I3I h32 /I33 

with hij being 3 by 3 matrices and H12 is given by (2.20). The result for the right boundary term is 

(0) + £^5^2(0) + ^7fc^3(0) ^ _ ^5fc^i(0) + ^4^2(0) + hkU^) 
k+{0) + ek+{0) + k+{0) ' fc^i(0) + efc2+(0) + fc3+3(0) ' 

zrk+jO) + ez.k+.jO) + z,k+{0) 

fc+(0) + efc+(0)+43(0) ^ ■ ^ 

The factor a in (2.23) is due to the normalization of H in (2.19) and e's appear in (2.26) to take into 
account the graded traces. Therefore the most general diagonal boundary terms can be written as 

b.t. = l{ln-h3)S! + l{ln+h3)iStr + l22li 

+ l{rn - r'33)Sh + ^ir'n + r'33)iSf,)^ + rasljv, (2.27) 

where /jj = la — I22 and r^^ = ru — r22 for i = 1,2,3. 

3 The coordinate Bethe Ansatz 

In this section results are presented for a open quantum spin chain of N atoms each with spin 1 described 
by the Hamiltonian (2.21) with the boundaries term (2.27): 

H =J2Hk,k+i+h.t. (3.1) 

k=l 

At each site, the spin variable may be +1, 0, —1, so that the Hilbert space of the spin chain is H^^^ = (g)^V 
where V = with basis {|+) , |0) , |— )}. The dimension of the Hilbert space is dimH^^^ = 3^. We can 
see that H commutes with the third component of the spin 

N 

[H,S^]=0, 5J = ^5fe^ (3.2) 

k=l 

This allows us to divide the Hilbert space of states into different sectors, each labelled by the eigenvalue 
of the number operator r = N — . We shall denote by Hn^'' the subspace of H^^^^ with r = n. We 
can see that dimH^^I = J2r=o dimif^"' with 

r-2j J \ j 



3=0 
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where [|] means the integer part of | and ^ 5 J denotes the binomial number. 
The action of -fffe,fe+i on two neighboring sites is read directly of (2.20) 

Hk,k+i = zi , Hk,k+i I — ) = zi I — ) , 

Hk,k+i\+0) = Z5|+0) + |0+), Ffc,fc+i |0+) = Z5 |0+) + 1+0) , 
Hk,k+i\0-) = Z5|0-) + h0), i?fe,fc+i|-0)=25|-0) + |0-), 

Hk,k+l\ + -) = ^7| + -) + 56|0 0) + ^3|-+), 

^ffe,fc+i|0 0) = e2:4|0 0)+ez6|+-) + ez6|-+), 

Hk,k+l\-+) = ^7|-+)+ -26|0 0)+^3|+-), (3.4) 

and the boundary terms (2.27) only see the sites 1 and A'' . Therefore, we can write 

b.t. |i, j) = £ij\i,...,j) , (3.5) 

where Sij ~ la + rjj, i,j = 1,2,3 , with the notation (+,0,—) = (1,2,3). la and rjj are given by 
(2.23) and (2.25), respectively. 

3.1 Sector r=0 

The sector Hq^^ contains only one state, the reference state, with all spin value equal to +1, ^'o = 
rife l+>^)> satisfying iJ^'o = Eq^q, with Eq = {N — l)zi + £\\. All other energies will be measured 
relative to this state. It means that we will seek eigenstates of H satisfying {H — £'o)^r — Sr'^r , in 
every sector r. 

3.2 Sector r=l 

In h[^\ the subspacc of states with all spin value equal to +1 except one with value 0. There are N 

AN) 

uj/ — -| — I — T"-] — r ■ ■ ■ -r / wiiiuii apaiL a uaaio ui iJ 
k I 

the form 

N 

*i=^a(fc)|fc[0]). (3.6) 

fe=i 

The unknown wavefunction a{k) determines the probability that the s'pin variable has the value at the 
site k. 

When H acts on |fc[0]) , it sees the reference configuration, except in the vicinity of k, and using (3.4) 
we obtain the eigenvalue equations 

{Si+2zi-Z5-Z5)a{k)=a{k-l)+a{k + l), {1< k < N) (3.7) 

At the boundaries, we get slightly different equations 

{£1 + £n - S21 + zi - Z5)a{l) = a(2), 

{£1 + £n - £12 + zi - -Z5)a{N) = a{N-l). (3.8) 

We now try as a solution 

a{k) = a{e)^'' - a(-6')^^ (3.9) 



states |A;[0])= + + +0 + + ••• + ) which span a basis of H{ ' . The Ansatz for the eigenstate is thus of 
fc / 
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where = e^^ , 6 being some particular momentum fixed by the boundary conditions. Substituting this 
in equation (3.7) we obtain the eigenvalue 

£^=-2zi + Z5 + Z5+C + r'^- (3.10) 

We want equations (3.7) to be valid for A; = 1 and k = N also, where a(0) and a{N + 1) are defined 
by (3.9). Combining (3.8) with (3.7) we get the end conditions 



o(0) = Aia(l), Ai = £21 - Sii + zl - Z5 = zi - Z5 - 

a{N + 1) = A2a{N), A2 = S12 - Sn + zl - Z5 = zi - Z5 - r'n. 
Compatibility between the end conditions (3.11) yields 

a{e) _^_2 Ai-^ _,_2N^2-^-' 



(3.11) 



= r^l=^=r^^=tr^, (3.12) 



a{-e) " Ai-^-i ^ A2-^ ' 

or 



i"=l¥^)(¥^}^ (3.13) 



Ai - e ; V A2 - e 

Therefore, the energy eigenvalue of H in the sector r = 1 is given by 

Ei = {N- 3)zi + hi +rii+Z5 + Z5+^ + (3-14) 
with ^ being solution of (3.13). 

3.3 Sector r=2 

In the Hilbert space H^^^ we have N states of the type \k[-]) = + + - + -\ h) and N{N - l)/2 

states of the type |fci[0], fc2[0]) 



+ + + +H h). We seek these eigenstates in the form 

fel fe2 / 

JV 

^ a(fci,fc2)|fci[0],fc2[0])+^6(fc)|fc[-])- (3.15) 

ki<k2 k=l 

Following Bethe [4], the wavefunction 0(^1,^2) can be parametrized using the superposition of plane 
waves (3.9) including the scattering of two pseudoparticles with momenta 9i and 02, (Cj = e*^^ , j = 1, 2): 

a{ki,k2) = ^ep{a{eue2)^^'^t-<02,e^)^^^^t} , (3.16) 
p 

where the sum extends over the negations of 61 and and ep is a sign factor (±1) that changes sign on 
negation. The paramctrization of b{k) is still midctcrmincd at this stage. 

Before wc try to parametrize b{k) let us consider the Schrodinger equation H'i>2 — E2'^2 ■ Prom the 
explicit form of H acting on two sites (3.4) we derive the following set of eigenvalue equations: 

• Equations for |fci[0]) and |fc2[0]) far in the bulk (1< fci < fc2 + l< A'') 

{£2+4:Zi-2z5-2z5)a{ki, k2) = a{ki- 1, fc2)+a(fci+ 1,^2)+ a(A;i, fc2-l)+a(fci, ^2 + 1). 

(3.17) 
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• Equations for \k[-]) in tlic bulk (1 < fc < iV) 

{£2 + 2zi - zr- zr)b{k) = Z3b{k - 1) + Z3b{k + 1) + ZQa{k -l,k) + z&a{k, k + 1). 

(3.18) 

• Equations for two |A;[0]) neighbors in the bulk {1 <k < N — 1) 

{£2 + 2>zi -z^-zr,- €Z4)a{k, k + 1) = a{k - 1, k + 1) + a{k, k + 2) + €Z6b{k)+ ezQb{k + 1). 

(3.19) 

In addition we have seven conditions to be satisfied at the free ends of the chain: 

• Five equations involving at least one state |A;[0]) at one of the ends 

{£2 + £n-£2i + 3zi-2z5- Z5)a{l, ^2) = a(2, ^2) + a(l, fe- 1) + a(l, ^2+ 1), 

(3.20) 

(^2+<?ii-^i2+ 3-21- Z!i-2-Z5)a{ki,N) = a{ki-l, N) + a{ki+ 1, N) + a{ki, N- 1), 

(3.21) 

{£2 + £11 - £22 + 221 -Z5- h)a{l, N) = a(2, N) + a(l, N - 1), (3.22) 
{£2 + £11 - £21 + 2zi -Z5- ez4)a{l, 2) = a(l, 3) + ez6b{l) + ezQb{2), (3.23) 
(£2+^11- £12 + 221- zr,-eZi)a{N- l,N) = a{N- 2, N) + ezeb{N- 1) + ezeb{N). 

(3.24) 

• Two equations with the state \k[—]) at one of the ends 

{£2 + £11 - £31 + zi- Z7)b{l) = Z3b{2) + Z6a{l, 2), (3.25) 

{£2 + £11 - £13 + zi- Z7)b{N) = Z3b{N - 1) + -zea{N - 1, N). (3.26) 
By simples substitution the Ansatz (3.16) solves the equations (3.17) provided 

£2 = -4zi + 2z5 + 2-Z5 + a + + 6 + ^2"'- (3-27) 

It immediately follows that the eigenvalues of H are a sum of single pseudoparticle energies. 

The parametrization of b{k) can now be determined in the following way: subtracting Eq.(3.19) from 
Eq.(3.17) for ki = k,k2 = k + 1, we get 

ezeb{k + 1) + sz6b{k) = a{k, k) + a{k + l,k + l) — {zi + SZ4 — Z5 — Z5)a{k, k + 1). 

(3.28) 

for which we can find b{k) in terms of a{ki, ^2)- 

Using (3.28) together with (3.17) we can see that (3.23) and (3.24) are readily satisfied. Now we 
extend the Ansatz (3.16) to fci = fc2 = A; in order to get a parametrization for the wavefunction b{k) : 

b{k) = Y.^lb{0u02)ili2 

p 

= b{eu92)ei^l- b{-e^,e2)^^H^ - &(^i -^2)^? ^2"' + ^-01-02)^^^2^ (3-29) 



which solves the meeting condition (3.28) provided 

A = Z5 + Z5 — zi — ez4, (3.30) 

together with b{—9i, 62), b{9i, —62) and b{—6i, —62) that can be obtained from (3.30) changing the signs 
of 01 and 62 ■ 

These relations tell us that the pscudoparticle of the type |fc[— ]) behaves under the action of H as 
the two pseudoparticles |fci[0]) and |fc2[0]) at the same site k and its parametrization follows as the plane 
waves of pseudoparticles \ki[0]) multiplied by the weight functions b{±9i,±92)- 

As a consequence of this identification we can see that the equation involving b{k) (3.18) becomes 
a meeting condition for two states |fc[0]). Using the S'-matrix language, from (3.18) we get the two- 
pseudopartcle phase shifts: 

where 

5(^2,^1) = (i+a6 + Ae2)[^3(i + eiC2)-(i+66)(a+6)+Aa6] 

+ £^2 (Ze + ^666) iz6 + 2666) I 

A = 2(zi - Z5 - Z5) + Z7 + Z7. (3.32) 

The seven remained s-functions for the phase shift equations (3.31) follow from (3.32) changing the signs 
of 9i and 92- 

At this point we still have to consider the equation (3.22) and four end conditions. We want equation 
(3.17) to be valid for ki = 1 and k2 = N also, where a(0, fc2) and a{ki,N + 1) are defined by (3.16). 
Combining (3.20) and (3.21) with (3.17) we get two end conditions 



(3.33) 



Aia(l, k) = a(0, k), Ai = zi - Z5 - £11 + £21, 
A2a{k, N) = a{k, N + 1), A2 = - - ^ii + ^12 
Substituting (3.16) in (3.33), we obtain the following relations 

^ r (^f^) (^) (3.34) 

which describe the change of signs of ^1 and 02 in a{6i,02) and the corresponding pair interchange 
relations 

t2Ar ^2 - 6 \ „//) /) X „f n a \ ^ ~ 



a(^2,-^i) = ^1 [ ^^ _^-i )ai(^2,0i), a{-02,0i)=l^^—^^j 0(02,01), 
ai-02,-0.) = 6- (^^) (^) ai02,0.). (3.35) 
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Combining these relations with the phase shift relations (3.31) we get the Bethe equations 



2^ _ n-AiCA f}_Ml\ {'JM)_\ f s{92,~er) \ 

2^ _ / i-Ai6 \ / i-A26 \ f s{e2,e,) \ f s{e„-0,) \ 

With these relations we have defined the behavior of the state |A:[0]) at the boundaries. Consequently the 
equation (3.22), where wo have one state |fc[0]) at each end, is also readily satisfied. 

Now we recall step by step our procedure described above to see that there is no more function or 
parameter to be determined, but we still have to solve the equations with the pseudoparticle |— ) at the 
boundaries. 

Similarly, we want equation (3.18) to be valid for k = 1 and k = N also, where 6(0) and b{N + 1) are 
defined by (3.29). Combining (3.25) and (3.26) with (3.18) we obtain two further end conditions 

A36(l) = Z36(0) + zea{0, 1), A3 = - ^7 - £11 + £31, 

(3.38) 

Aih{N) = Z3b{N + 1) + zea{N, N + 1), A4 = zi-Z7- Sn + S13. 

Substituting the already fixed relations for a{ki, ^2) and b{k) and using the Bethe equations we can see 
that these end conditions are also satisfied. It means that there is no additional end condition due to the 
presence of the state \k[—]) at the boundaries. Instead of surprising, this result is in agreement with our 
paranictrization of the wavefunction b{k), where the dynamics of the pseudoparticle \k[—]) is understood 
as the dynamics of two pseudoparticles |A;[0]). Therefore, the energy eigenvalue for the sector r = 2 is 

2 

E2 = {N- l)zi + In +rii+J2 ("^^i + ^5 + h + + , (3.39) 



with 



f2JV 



1-AigA / i-A2gA f sie,,0k) \ f si9k,-e,) 



j = 1,2. (3.40) 
where the functions s{9j,9k) are given by (3.32). 

3.4 General sector 

The above results can now be generalized to arbitrary values of r. In a generic sector r we build 
eigenstates of H as direct products of A^o states |/c[0]) and N- states |A;[— ]), such that r = Nq + 2N- . 
These eigenstates are obtained by superposition of terms of the form 

\4>r) = |0) X + I-) X \^r-2) , (3.41) 

with l^o) = 1) |0i) = |0). For instance, in the sector r = 3 the eigenstate of H has the form 
*3 = J2 aikuk2,k3)\ki[0],k2[0],k3[0]) 

fel <fe2 <fe3 

+ ^ {bl{ki,k2)\ki[-],k2[0])+b2{ki,k2)M0],k2[-])}. (3.42) 

fel <fe2 
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The Ansatz for the wavefunction of the term with iVo states |fc[0]) becomes 

a(fci,fc2,... ,M=E^f«(^l'^2,... ,^r)Cl'^C2'---^r'-> (3-43) 
P 

where the sum extends over all permutations and negations of ^i, ^2, • • • ,dr and Ep changes sign at each 
such mutation. 

The Ansatz for the wavefunction of terms with N- states \k[—\) follows from (3.43) as sum over 
negations of the terms with 2N- states |A;[0]) at the same site. For instance, in the sector r = 3 



13(^2, 03)^1^ e2^j^l^ (3.44) 

with similar equation for &2(^ij ^2). In that way we always have a far pscudoparticle |fc[0]) as a viewer. 
We also have verified in this sector that the meeting of |fc[0]) with \k[—]) can be versed as a meeting of 
three |fc[0]) whose phase shift factorizes in a product of two-pseudoparticle phase shifts. 
The corresponding energy eigenvalue is a sum of single one-particle energies 

r 

Er^{N- l)zi + hi + rii + E (-2zi + Z5 + Z5 + + ij^) , (3.45) 
where are solutions of the Bethe equations 

j = l,...,r (3.46) 

s{ej,ek) = (1 + + A?,) [^3(1 + e -eD - (1 + ^j^km + ^k) + A^^-a] 

+e^j {Z6 + ZQ^jik) {zq + ze^j^k) , (3.47) 

A = 2{zi - 25 - Z5) + Z7 + Zy, 

Ai = 2:1-25-^11+^22, A2 = zi - Z5 - ?'ii + ?'22. (3.48) 

Having now built a common ground for all open spin-1 Hamiltonians associated with the 19-vertex 
models, we may proceed to find explicitly their spectra. We will do that in the next sections. 



with 



and 



4 The Zamolodchikov-Fateev model 

The simplest three-states vertex model is the ZF 19-vertex [15] or the A^^^ model the spin-1 representation 
[20] and can be constructed from the six-vertex model using the fusion procedure. The 7?.-matrix which 
satisfies the YB equation (2.1) has the form (2.18) with 

xi{u) = sinh(w + 77) sinh(u + 277), X2 (w) = sinh it sinh(M + 77) , 

2^3 (u) = sinh u sinh(M — 77) , X4 (u) = sinh u sinh(u + r]) + sinh 77 sinh 277, 

7/5 (u) = X5{u) = sinh(7i + ry) sinh27y, yeiu) = xe{u) = sinh u sinh 2r/, 

y7{u) = a;7('u) = sinh7/sinh2r/. (4.1) 
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This 7?,-inatrix is regular and unitary, with f{u) = xi{u)xi{—u), P- and T-symmetric and crossing- 
symmetric with M = 1 and p = r}. The most general diagonal solution for K~{u) has been obtained in 
Ref. [19] and is given by 




i^-(w,/3n)=| 1 1, (4.2) 
with 

riiuh » + 2 cosh u, .-/>,_ siuh((( + ;/) - 2 cosh^;; + (7) 

''^"^""/3iisinhu-2coshu' 33(«) - -^^^ - J?) + 2 cosh(« - j?) ' ^ ^ ^ 

where /3ii is the free parameter. By the automorphism (2.7) the solution for K^{u) follows 

K+{u,ari)=K-{-u-p,a^^)=\ 1 |, (4.4) 



with 



, + / X Q!ii sinh(u + r/) — 2 cosh(zt + 77) , + . x an sinhu + 2 coshw 

fcii(w) = , „^ , o„„„u/„. , „^ ' ^33(«) = 



a\i sinh('u + r/) + 2 cosh(u + r?) ' an sinh(u + 2r/) — 2 cosh(u + 2vi) ' 



(4.5) 



where an is another free parameter. 

We recall section 2 to derive the corresponding quantum open spin chain Hamiltonian. It is the 
quantum open spin chain for the spin-1 XXZ model. 

N-l 

H=Y,Hk,k+i+h.t. (4.6) 
fc=i 

where the bulk Hamiltonian is given by (2.21) with the weights 

e = 1, a = sinh2ry, 2:1 = 0, 23 = —1, 2:4 = — 2cosh2r7, 
Z5 = Z5 = — cosh 277, Z6 = ze = 2cosh?7, zr = Z7 = —1 — 2cosh2r7. (4.7) 

and the boundary terms given by (2.27). For the left boundary (2.23) we get 

1 ^ . , „ , „ , /3n cosh 71 — 2sinh?7 . , „ , ^ 
In = -I3nsmh2r], ^22 = 0, ^33 = ? ■ . o r^sinh277, 4.8 

2 /Jii smhjy — 2cosh77 

and 

ail coshr; — 2 sinhr? . , „ 1 . , „ 

r-ii - r22 = r-r 7, r- smh2r?, r^s - r-22 = ;raii smh2?7, 

an smh ?7 — 2 cosh ?7 2 

lsinh477 / [ail sinhjy + 2cosh?7]^ — 4[1 + 2cosh2r7]\ (4 9) 

4sinh3ry\ an sinhry — 2coshr7 /' 

for the right boundary (2.26). 

Next, we can use the coordinate BA method as described in section 3 to find the energy eigenvalue 
(3.45) and the BA equations (3.46). Here we just list the results. The energy spectrum of the Hamiltonian 
(4.6) for a generic sector r is given by 

„ /I an cosliT) — 2sinh77\ -sr^ , , ^ ^ is , 

E. = sinh2, ( -Ai + ] + + E(-2cosh2, + + ^7^), (4.10) 
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with $,j = e ' satisfying the Bcthc equations 



c2N 



l,2,...,n (4.11) 



where 



A » ail coshr? — 2sinh?7\ 

Ai = sinh277 coth27? - -jSn ] , A2 = sinh2n coth2n . ' — — ^ . 

\ 2 J \ ail sinhTj — 2cosh77/ 



For the ZF model the two-particle phase shift is given by 

s{ej,9k) _ 1 + + g + gj-g - (a + 2)^^- 



siOk, Oj) l + ^j+^k+ O^fc - (A + 2)g ' 



(4.12) 



(4.13) 



with the s-functions given by (3.47) and A = 2 cosh 277. 

This energy spectrum was already obtained by Mezincescu at al. [19] through a generalization of the 
quantum inverse scattering method developed by Sklyanin [12], the so-called fusion procedure [20]. This 
fusion procedure was also used by Yung and Batchelor [21] to solve the ZF vertex-model with inhomo- 
geneities. 

For a particular choice of boundary terms, the ZF spin chain has the quantum group symmetry i.e., 
if we choose ^ 00 {Pu = 2coth^_ and an = 2coth^+), then the spin chain Hamiltonian (4.1) has 
i7q(su(2))-invariance [19]. 

5 The Izergin-Korepin model 

The solution of the YB equation corresponding to A2 in the fundamental representation was found by 
Izergin and Korepin [16]. The 7^-matrix has the form (2.18) with non-zero entries 

Xi{u) = sinh(u — 677) + sinhry, a;2(u) = sinh(M — 3?]) + sinh3r/, 

2^3 (m) = sinh(u — 77) + sinhry, a;4(u) = sinh(M — 37^) — sinhSry + sinhSry + sinhr; 

X5{u) = -2e~"/^sinh 277Cosh(^ - 3??), y5(w) = -2e"/^ sinh 2?] cosh(^ - 3??) 

xe{u) = 2c-"/2+2''sinh 277 sinh(^), ?/6(«.) = -2e"/2-2'; sinh277sinh(|) 
xt{u) = — 2e^"^^'' sinhr7sinh2?7 — e"** sinh47y, 

Uriu) = 2e"~^'' sinh77sinh 277 — e'' sinh477. (5-1) 

This 7?.-matrix is regular and unitary, with f{u) = Xi{u)xi{—u). It is PT-symmetric and crossing- 
symmetric, with p= — 67] — iiT and 

M = 1 . (5.2) 

Diagonal solutions for K~{u) have been obtained in [22]. It turns out that there are three solutions 
without free parameters, being K~{u) = 1, K^{u) = F+ and K~{u) = F~, with 

V(±)(w) 

(5.3) 

e"/(±)(u) 
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where we have defined 

.(±) ^ ^ ^ cosh(u/2 -3T])±i sinh(u/2) 

^ ^ cosh(u/2 + 377)Tisinh(u/2)' ^ " ' 

By the automorphism (2.7), three solutions K~^{u) follow as K~^{u) = M, K^{u) = G+ and K^{u) = G~, 
with 

/ e"-4''5(±)(u) \ 
= 1 , (5.5) 

where we have defined 

J±) (,A ^ cosh(M/2 -3T])±i sinh(M/2) 

^ ^ ' cosh(M/2-3?7)T«sinh(u/2-6r/)' ^ ' 

The corresponding quantum open spin chain Hamiltonians is also written as in (4.6), where the bulk 
term is given by (2.21) with 

„ . , „ „ cosh 7? „ sinh 77 sinh 4r7 

e = 1, a = -2smh277, = 0, 03 = r^, Zi = -2 '— '- 

cosh 3?7 cosh 877 

_2„ . 2„ 2r7sinh277 . _2„sinh277 

Z5 = -e 2:5 = -6 ", -26 = 6^'' — — — , ze = -e — — — 

cosh 377 cosh 37; 

= _^^(e-4'7 + 2sinh27?), ^7 = - (e^" - 2 sinh 27?) . (5.7) 
cosh 377 ^ ' cosh 377 ^ ' 

To derive the boundary term (2.27), we will only consider three types of boundary solutions, one for each 
pair {K-{u),K+{u)) defined by the automorphism (2.7): (1,M), {F+,G+) and (F-,G-). 

The (1,M) case: For K-{u) = 1 the left boundary (2.23) vanishes ( hi = I22 = hs = ), 
and the right boundary (2.24) is proportional to the identity. To see this we can substitute (5.7) and 
K+{0) = diag(e2'', l,e-2')) in (2.26) to get 

„cosh47?sinh277 _ 

rii = r22 = r33 = -2 r^- '-. 5.8 

sinh 6?7 

Therefore, the corresponding open chain Hamiltonian is 

7V-1 



cosh 4r; sinh 21] ^ 
sinh 677 

The coordinate BA gives us the corresponding energy spectrum (3.45). For a given sector r it is given 



H^J: H,,,^i - 2 ^"---;-'/ — -/ i, (5.9) 
smn 677 

fc=i 



by 

sinh 677 



^ ^ cosh 477 sinh 277 / „ , „ ^ ^ i\ 

= -2 Aar. + E (-2 cosh2,7 + + ^-1) , (5.10) 



where = e'^^ are solutions of the Bethe equations 

c2N _ TT f s{0j,9k)\ f s{9k,-0j] 



r = n 



j = l,...,r (5.11) 
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The two-particle phase shift for the IK model is given by 



s{ek, 0j) V 1 + -A^kJ \ i + Cj-a - - a - (A - 2)^. 



(5.12) 



where the s-functions are given by (3.47) and A = 2 cosh 2??. 

It was noted by Mezincescu and Nepomechie [22] , that this open spin chain Hamiltonian is quantum 
group invariant. Moreover, the corresponding transfer matrix has been diagonalized by the analytic BA in 
[23], using the Uq{su{2)) invariance of (5.9). This quantum group invariance was also used by Yung and 
Batchelor [21] to determine properties of the transfer matrix eigenvalues with inhomogeneities essential 
to apply the analytical BA. 

The {F~^, G^) and {F~,G~) cases: These cases can be treated simultaneously. The matrix elements 
of the left boundary term are given by 

and for the right boundary term we have 

rW-rg) = -sinh2,(i^), - 4? = «inh2, 

(±) sinh4ry /coshTr? ± 4isinh37ysinh?7sinh2?7\ 

sinh 6rj \ cosh 3ri±i sinh 2ri J 

In these cases the corresponding open chain Hamiltonians are not J7g(sw(2))-invariant [23]. Nevertheless, 

it has recently been argued by Nepomechie [24] that the transfer matrices corresponding to these solutions 
also have the t/g(o(3)) symmetry, but with a nonstandard coproduct. They can be written in the following 
form 

From the coordinate BA we have find their energy spectra (3.45): 

= 2sinh2^ ( ''eosSr ) + "''^ + ^"^ + + ^^'^ ' ^^'^^^ 



where = e^^ are solutions of the Bethe equations 



) [ Al-) A^^ - ) .ii, Uom) U-Oj,0, 

3 = l,...,r (5.17) 



where 



A'«=e»--sinh2,(^), A<« = + stoh 2, ( ^) . (5.18) 

and the two-particle phase shift is still given by (5.12). These cases were also considered in Ref.[21] 
through the analj^tical BA with inhomogeneities. 
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Finally, we note that is interesting to reformulate the Bolthzmann weights of the IK model by the 
following transformation 

n{u, 7]) ^ n {u, 7]) = ^n{2u, -v-i^)- (5.i9) 

This TZ matrix differs from the one given in [25] by a gauge transformation. It is regular and unitary, 
with / (u) = Xi{u)xi{—u), PT-symmctric and crossing-unitarity with M = diag(— e"^**, 1, — e^'') and 
p = 3r]. After the gauge transformation TZi2{u) = ViTZi2{u)V^^ with V = diag(e~", l,e"), can see that 
M = diag(— e^'', 1, — e"*'') and p = p . In this case the solution G+) can be written as 

y'- =diag(l,-!!^^^|!i^,l), G"- = -diag(e^^!!^^^|^,e-^). (5.20) 

This solution was used by Fan in [26] to find the spectrum of the corresponding transfer matrix using 
the algebraic BA for one and two-particle excited states. 



6 The s/(2|l)-model 

The solution of the graded YB equation corresponding to s^(2|l) in the fundamental representation has 
the form (2.18) with non-zero entries [27, 17]: 

Xi{u) = cosh{u + r]) smh{u + 2r]) , a;2('«) = sinhucosh(M + r/), 

xaiu) = sinhucosh(u — 77), X4{u) = sinh'Ucosh(u + r]) — sinh277C0sh77. 

2/5 («) = X5{u) = sinh 277 cosh(u + 77), ye{u) = xe{u) = sinh 277 sinh «, 

yriu) = a;7(u) = sinh277cosh77. (6.1) 

This 7^-matrix is regular and unitary, with f{u) = xi(m).xi(— w,), P- and T-symmctric and crossing- 
symmetric with M = 1 and p = rj. The graded version of the crossing-unitarity relation (2.4) is satisfied 
with / (u) = x-i{u + i^)xi{—u — i^). 

The most general diagonal solution for K~{u) has been presented in Ref. [28] and it is given by 




K-{u,Pii)=\ 1 |, (6.2) 
with 

^_ , , /3ii sinhu -I- 2coshu ^_ /Jn cosh(u -I- 77) — 2 sinh(u, + 77) (Q3) 

I3n sinh u - 2 cosh u ' 33V"; - cosh(w - 77) + 2 sinh(M - 77) ' 

where /3ii is the free parameter. Due to the automorphism (2.7) the solution for K~^{u) is given by 
K~{-u - p, iaii) i.e. 

( kUu) \ 
K+{u,l3r^)=\ 1 , (6.4) 

V ^sW / 

where 

^+ / \ '^11 cosh(7i + 7/) — 2 sinh(?i + 77) ( ) '^^^ smh u + 2 cosh u 

"^"'^ ~ aiiC0sh(77-F7/)-F2sinh(77-^77)' 33 W - q,^^ sinh(77 277) - 2 cosh(7x 2?/) ' 
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and ail is another free parameter. 

The weights for the corresponding bulk Hamiltonian (2.21) are given by 

e = —1, Q: = sinh2?7, = 0, Zs = 1, ^4 = 2cosh2?7, 
Z5 = Z5 = — cosh2r/, zq = zq = 2sinh77, zr = z^ = 1 — 2cosh277 (6.6) 

The left boundary terms of b.t. (2.27) are given by 

; • ; n ; sinhr; - 2 coshr? . 

hi = 7r/3ii smh27?, I22 = 0, ^33 = . smh2j7 (6.7) 

2 fill cosh T] — 2 sinh 77 

and for the right boundary we have 

ail sinhn — 2 coshr? . , „ 1 . , „ 

ni - r22 = r „ . , smh2r?, r^s - r22 = ^an smh2?7, 

ail cosh T] — 2 smh 77 2 

1 sinh 477 /(ail cosh + 2 sinh 77)^ + 4(1 - 2 cosh 27?) \ 

T22 ~ I ) • (^•^) 

4 cosh 377 \ ail cosh 77 — 2 sinh 77 / 

Now, using the coordinate BA we find the energy spectrum (3.45) for the s/(2|l) open chain Hamil- 
tonian; 

^ . , „ /l^ ail sinh?7 — 2 cosh?7\ v^/ ^ , „ ^ , 1n „x 

= smh 277 [-,Pu + )+r2.+ E(-2cosh277 + + ^7^) (6-9) 



with = e'^' satisfying the Bcthc equations 



j = 1,2,. ..,71 (6.10) 

where 

1^^ A •uo/^ii,o„ ail sinh 7/ — 2 cosh 77^ 

(6.11) 



.,„ / \ . .,„ / ,„ ail sum 7/— ^COBU7/ \ 

Ai = sinh 277 coth 2r]- -Pu] , A2 = sinh 277 coth 27? —!■ — ^ 

\ 2 / \ ail cosh 77 — 2 smh 77/ 

and the two-body phase shift for the sZ(2|l) model is given by 

s{ej,9k) ^ l-g, -a+gj-a-(A-2)$, 

s{ek,ej) i-^j-a + eia-(A-2)a 

where A = 2 cosh 277. 



(6.12) 



7 The osp( 1 1 2)-model 

The trigonometric solution of the graded YB equation corresponding to osp{l\2) in the fundamental 
representation has the form (2.18) with non-zero entries [17]: 

xiiu) = sinh(u + 2r/) sinh(M + 877), X2(w) = sinh ?y, sinh(M + 877) 

xsiu) = sinh u sinh (u + 77), 2:4(71) = sinh Msinh(M + 877) — sinh 2?7 sinh 877 

3^5 (u) = e~" sinh 277 sinh(u + 877), 775 (it) = e" sinh 27/ sinh(w + 877) 

xq{u) = —e~"~^'' sinh 277 sinh 7x, j/e (w) = e""'"^'' sinh 277 sinh u 

Xt{u) = e~" sinh 277 (sinh('U + 877) + e"'' sinh u) 

yriu) = e" sinh 277 (sinh(u + 377) +e'' sinh m) (7.1) 
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This 7?,-niatrix is regular and unitary, with f{u) = xi{u)xi{—u). It is PT-symmetric and crossing- 
symmetric, with p = 3r] and 

M = 1 . (7.2) 

Diagonal solutions for K~{u) have been obtained in [29]. It turns out that there are three solutions 
without free parameters, being K^{u) = 1, K^{u) = and K~{u) = F~, with 

F± = 1 , (7.3) 

where we have defined 

fi+) = sinh(^ + 3r//2) fi-)r, cosh(u + 377/2) 

^ ' sinh(u-3r//2)' ^ ' cosh(u - 3r//2) ' ^ ' 

By the automorphism (2.7), three solutions K~^{u) follow as K'^{u) = M, K'^{u) = G+ and K~^{u) = G~ , 
with 

/ ^e2"+%(±)(«) \ 
G± = 1 , (7.5) 

\ ^g-2«-4r,^(±)(^) j 

where we have defined 

„(+) _ sinh(K + 377/2) (_) _ cosh(» + 377/2) 

^ ^""^ ~ sinh(w + 9r//2) ■ ^ ~ cosh(u + 97?/2) ' ^^-^"^ 

The corresponding quantum open spin chain Hamiltonians is also written as in (4.6), where the bulk 
term is given by (2.21) with 

. , „ „ sinh77 „ cosh 77 sinh 477 

e = -1, a = smh2?7, Zi=Q, ^3 = . , o ' ^4 = 2 —- 

smh 3?7 sinh 3?7 

2 . _2« _2„sinh2r/ . 2«sinh2?7 

Z5 = -e Z5 = -e ', ZQ = -e ' . , „ , 25 = 6 . , „ , 

smh 3 7/ smh 3?? 

2„ _„ sinh 277 . _2„ „ sinh 277 . . 

To derive the boundary term (2.27), we will only consider three types of boundary solutions, one for each 
pair {K-{u),K+{u)) defined by the automorphism (2.7): (1,M), (F+,G+) and (F-,G'-). 

The (1,M) case: For K~{u) = 1 the left boundary (2.23) vanishes ( In = I22 = ^33 = ), and the 
right boundary (2.24) is proportional to the identity, for which quantum-algebra invariance is achieved 
[23]. To see this we can substitute (7.7) and ^^+(0) = diag(e-2'', 1,6^'') in (2.26) to get 

cosh 477 sinh 277 .7o^ 

rii = r22 = 7-33 = 2 —— . (7.8) 

smh 677 

Therefore, the corresponding open chain Hamiltonian is 

R V^'h- I ^ cosh 47/ sinh 277 

H=}_^ + 2 liv (7.9) 
fe=i ' 
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The coordinate BA gives us the corresponding energy spectrum (3.45). For a given sector r it is given 

by 

^ „ cosh 471 sinh 2r? / i\ 

= ^ sinh6^ + ^ + + ^^-^^^ 

where = e'^^ are solutions of the Bethe equations 

j = l,...,r (7.11) 



The two-particle phase shift for the osp{l\2) model is given by 

s{0j,ek) ( 1 + ^j^k - A^A / 1 + ^jCfc + + Cfc - (A + 2)a 



s(^fc, 0,) V 1 + ^jik - Aa ; V 1 + ^j^fc + ^3 + a - (A + 2)^,- 

where the s-functions are given by (3.47) and A = 2cosh2r/. 
The (i^+,G+) case: In this case the boundary terms are 



(7.12) 



C-Sn/Z g3»7/2 

= • y^to /o^ sinh2?7, ^22 = 0, = . , , , sinh2?7 (7.13) 
smh(3ry/2) smh(3ry/2) 



and 



= sinh(3^/2)^^"^'^' "^^-"^^=sinh(3^/2)^^"^''^ 
The energy eigenvalues are 

r 

Er = 2sin2r?coth(3r//2) + r22 + ^ (-2cosh277 + + ^7^) (7.15) 

with the Bethe equations 

2JV ^ (^ 1-Ai0 \ / l-A2gA /£(£^\ /£(^fc^ 

V Ai - A A2 - e,- ; . 0,)J \s{-0j,0k) 

j = l,...,r (7.16) 

where the phase shift is given by (7.12) andAi = e-^") - ^^^t^^^^''^^'^) sinh2ry, A2 = e^'' - ^i^f^g^Jy^) sinh277 
The (F~, G~) case: In this case the boundary terms are 



e-Srj/S g3r,/2 

sinh277, ^22 =0, ^33 = sinhr? 



cosh(3r//2) " ' cosh(3r;/2) 

g3r)/2 g-3r,/2 

r-11 - r22 = — , /o^ ^^^^ ^'7' ''33 - r-22 = , /o\ ^^^^^ ^J? 

cosh(377/2) cosh(3?y/2) 
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The corresponding energy eigenvalues are given by 



Er = 2sin2r?tanh(3r?/2) + r22 + ^ (-2cosh2r/ + + ^r^) (7.18) 
The Bethe equations are 

j = l,...,r (7.19) 
with the phase shift (7.12) and 

Ai=e-^'> + — , A2=e2'' — — -sinh?? (7.20) 

cosh(377/2)' cosh(3?7/2) ' ^ ' 



8 From non graded to graded solutions 

Beside the 7?.-matrix we also have considered the i?-matrix, which satisfies 

Rl2{u)R23{u + v)Ri2{v) = R23{v)Ri2{u + v)R23{u). (8.1) 

Because only R12 and R23 arc involved, this equation written in components looks the same as in the 
non graded case. Moreover, the matrix TZ = PR satisfies the usual YB equation (2.1) where P is the non 
graded permutation matrix. When the graded permutation matrix V is used, then TZ = PR satisfies the 
graded version of the YB equation. 

Multiplying the 7^-matrix for 19- vertex models (2.18) by the diagonal matrix 11 = PP = PP wc will 
get graded 7?.-matrices starting from non graded 7?.-matrices and vice- versa. The new 7?.-matrix TZ = UTZ, 
still has the form (2.18) but with the change of sign of the fifth row due to the grading BFB. The bulk 
Hamiltonian has the form (2.20) but interchanging the role of the sign e. Now e = — 1 for non-graded 
models and e = 1 for graded models. 

Let us use this interchange property with the YB solution of the IK model. First we recall the 
transformation (5.19) 

TZ'{u,r]) = ^n{2u, -n - z|) ^ <fe+i(r,) = i?fe,fe+i(-ry - i|). (8.2) 

The matrix TZiKg{u,ri) = UTZ is a solution of the graded version of the YB equation (2.1) and the 
corresponding vertex model can be named as the graded version of the IK model. 

Using the symmetries of the YB solutions for 19- vertex models: X2{u) — > ±X2(m) and X(i{u) ±xe{u) 
with ye{u) TUeiu), we can see that this model has the same Boltzmann weights of the osp( l|2)-model, 
except for the presence of the factor ±i in xq{u) and in 2/6 (w). However, this identification is 
not so trivial due to the change the signs of the fifth row of TZ (BFB grading). Nevertheless, by direct 
computation we have verified that both models have the same reflection ii'-matrices. It means that 
'R-iKg{u,rf) and the TZ{u,rf) of the osp(l|2) share the same symmetries. Consequently, both open chain 
Hamiltonians have the same boundary terms. Moreover, from the definition (3.32) we can see that 
phase shift equations (3.31) are invariant under the replacement 2:5 iiz^ with zq T^Zq. Thus, the 
coordinate BA previously described, yields the same spectrum for both models. In words, the open spin 
chain Hamiltonians associated with the graded IK model have the osp{l\2) - invariance. 
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This situation is also present in the graded version of the ZF model. In order to see that we have to 
reformulate conveniently the Boltzmann weights of the ZF model by the following transformation 

n{u,7])^n{u,7]) = lTl{u,r,-il). (8.3) 

The graded version of the ZF model is defined by the following 7?.-matrix 

nzFg{u,v)=Tm'{u,ri) (8.4) 

Using again the symmetries of the 19-vertex model we can sec, up to a possible canonical transformation: 
xe — > Xg{u) =±ix(i{u), the non-zero entries of TZzFg{u,ri) are identified with the Boltzmann weights of 
the s^(2|l) model (6.1). We also find that both models have the same ii'-matrices and their coordinate 
Bethe ansatze yield a common spectrum. 

We have verified that the inverse situation is also true. The non-graded versions of the graded 19-vertex 
models are in correspondence with the 19-vertex models of Izergin-Korepin and Zamolodchikov-Fateev. 

During the preparation of this paper we learned that the connection between Izergin-Korepin and 
osp{l\2) models has recently been discussed in Saleur and Wehefritz-Kaufmann [30], where also earlier 
references are given. 

9 Conclusion 

In the first part of this paper we have applied the coordinate BA to find the spectra of open spin-1 chain 
Hamiltonians associated with four 19-vertex models, including two graded models. This procedure was 
carried out for boundaries derived from diagonal solutions of the reflection equations. 

We believe that the method here presented could also be applied for Hamiltonians associated with 
higher states vertex- models. For instance, in the quantum spin chain s = 3/2 XXZ model we have four 
states: |fc[3/2]), \k[l/2]), |A;[-l/2]) and |fc[-3/2]) . It means that the state |A:(l/2)) can be parametrized 
by plane wave and the states 1/2]) and 3/2]) as two and three states ]A;[l/2]) at the same site , 
respectively, multiplied by some weight functions. 

These weight functions are responsible by the factorized form of the two-body phase shifts of the IK 
model (5.12) and the osp(l\2) model (7.12). In the ZF model, as well as in the s^(2|l) model, we do 
not have a factored form for the two-pseudoparticle phase shift because their weight functions (3.30) are 
constant. It means that the state ]) behaves exactly as two states \k[0]) at the same site. This is in 
agreement with the fact that the ZF model can be constructed by a fusion procedure of two six- vertex 
models. 

There are several issues left for future works. A natural extension of this work is to consider the 
algebraic version for the BA [26]. Independently, it is interesting to analyse the Bethe Ansatz equations 
to derive ground state properties, low-lying excitations and the thermodinamic limit. 
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